Early stages of magnetization relaxation in superconductors 
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Magnetic flux dynamics in type-II superconductors is studied within the model of a viscous 
nonlinear difi'usion of vortices for various sample geometries. We find that time dependence of 
magnetic moment relaxation after the field is switched off can be accurately approximated by m{t) oc 
1 — -s/t/f in the narrow initial time interval and by m{t) oc (l + t/r)~^ at later times before the fiux 
creep sets in. The characteristic times f and r are proportional to the viscous drag coefficient 77. 
Quantitative agreement with available experimental data is obtained for both conventional and high- 
temperature superconductors with 77 exceeding by many orders of magnitude the Bardeen-Stephen 
coefficient for free vortices. Huge enhancement of the drag, as well as its exponential temperature 
dependence, indicates a strong influence of pinning centers on the flux diffusion. Notwithstanding 
the complexity of the vortex motion in the presence of pinning and thermal agitation, we argue that 
the initial relaxation of magnetization can still be considered as a viscous flux flow with an effective 
drag coefficient. 

PACS numbers: 74.25.-q, 74.25. Wx 



Magnetic flux penetrates a type-II superconductor in 
the form of discrete quantized vortices. Vortex struc- 
tures in conventional and high-temperature superconduc- 
tors display remarkable complexity both in equilibrium i.!^ 
and dynamic regimes i^"— Relaxation of the magnetic mo- 
ment of superconductors is achieved through initial vis- 
cous flux flo'W'i^"— and slow, logarithmic in time, ther- 
mally activated creep.-i§"-'^ Thermally-assisted hopping 
of vortices and vortex bundles between local minima in 
the random pinning potential is characteristic of both 
the creep and the flux flow under a driving force. In 
the latter, the hopping gives rise to the viscous drag 
coefficient 77 oc e'^^'^^, where U is the effective acti- 
vation energy and T is the temperature j2i A free flux 
flow regime can be realized at microwave frequencies 
(10 — lOOGHz) when the effect of the pinning is negli- 
gible. Measurements of surface impedance give viscous 
drag coefficients 770 ^ 10^^ — 10"^ Ns/m^ at low tempera- 
tures for all superconductors, e.g., conventional NbSe2r^ 
cuprates YBCO, BSCO^i^ and pnictide LiFeAsJ^ The 
order of magnitude is in accordance with the Bardeen- 
Stephen result for the viscous drag, 770 = ^oHc2/ PnC^ , 
caused by dissipation in the vortex core ($0 = hc/2e is 
the flux quantum, p„ is the normal-state resistivity, and 
is the upper critical field) 4. 

In this paper we study early stages of the fiux dy- 
namics after switching off the external magnetic field. 
We use a simple hydrodynamic approach: The local 
force the vortex experiences due to interaction with other 
vortices, the surface, and the local quenched disorder 
(pinning centers) is described by an effective viscosity 
'7 ^ ?7o.The same approach successfully describes the vor- 
tex creep, if supplemented by a phenomenological model 
of current-dependent or time-dependent activation en- 
ergy, U = Ucln{jc/j) OT U = kT\n{t/to), where jc is 
the critical current and Iq is the characteristic time scale 
for fiux creep J^^— 



We consider a model of massless vortex motion where 
the driving Lorentz force equals the viscous drag (1/c) J x 
^Q — r]v = 0. Here, J is the current density, v is the vor- 
tex velocity, and 77 is a viscous drag coefficient. For mag- 
netic induction B — 7i$o related to the vortex density n, 
the force balance equation reads (l/c)JxB — 77|i3|v/<I>o = 
0, with J = (c/47r)V x B. Taking into account the con- 
tinuity equation dB/dt + V • {Bv) = 0, the dynamics of 
the magnetic flux in a superconductor is described by the 
well-known nonlinear diffusion equation^Sr— 



dB $0 , . 

dt iTTT] ^' ' ' 



(1) 



We have solved Eq. ([T]) for three sample geometries: a 
slab, a square-shaped plate, and a disk (see Fig. [J). We 
assume the sample thickness along the field is sufficiently 
large and neglect stray fields on the top and bottom of 
the sample. Magnetic induction B(r,t) is directed along 
the sample symmetry axis z and satisfies the following 
initial and boundary conditions: (i) B is uniform within 
the sample at t = 0, B{x,y;t = 0) = Bq, and (ii) B 
vanishes at the sample edges for t > 0. 

For a long superconducting slab of width L, Eq. ([T|) 
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FIG. 1. Vortex dynamics is studied for three sample geome- 
tries: (a) a slab, (b) a square-shaped plate, and (c) a disk. 
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FIG. 2. Magnetic induction B{x) across one half of a super- 
conducting slab of width L is shown for times t/ro — 0.01, 
0.1, 0.2,. . . , 3.2 (top to bottom). 



reads 



dt 47r?7 dxV '9a;/' 



(2) 



where B{x,t) = at a: = ±L/2 for t > 0. We can seek 
the solution in the form 

oo 

Bix,t) = J2Bkit)sm[kTr{x/L + l/2)], (3) 
fe=i 

with functions Bk{t) to be determined from Eq. 1^ and 
B{x,t = 0) = Bq. This gives the following set of differ- 
ential equations: 



dBkjt) 
dt 



1 



BqTq 



B,{t)Fk{i,j)B,{t), 



(4) 



with the initial conditions Bk{0) = {2Bo/Trk) [1 - {-!)''] 
(fc = 1, 2, . . .). Here, the coefficients Fk{i,j) are given by 



(5) 



for |i±j| ^ k and i+j+A: odd, and Fk{i,j) = otherwise. 
The characteristic time constant is 



To 



$0^0 ■ 



(6) 



Equations ^ are solved by truncating the system at suf- 
ficiently large k [k ^ 40). The induction B{x,t) for the 
slab is shown in Fig. [2] at various times t/ro — 0.01, 0.1, 
0.2, 3.2. We observe that the flux flow near the sample 
edges in the initial time interval is very fast, reaching the 
center of the slab {x = 0) at time < ^ 0.1 tq after switch- 
ing off the field. This regime is followed by a slower flux 
flow taking place in the bulk of the sample. 

The spatial dependence of the magnetic induction is 
in accordance with the previous results for the flux flow 



FIG. 3. Numerical solutions for magnetic moment relaxation 
for the slab (x), disk (o), and square (□) geometries, com- 
pared to the analytic approximation, Eq. ((9| (solid line). In- 
set: A better fit for the initial time interval t <^ t to the 
analytic expression, Eq. ((7|. In this case, only magnetization 
at the edges is affected by the flux flow. 



regime with constant activation energy^S^ In the pres- 
ence of flux creep, which may take place in the cen- 
ter of the slab for t ^ O.Itq, or at large times t ^ tq 
when remanent magnetization is small, a phenomenolog- 
ical model of current and fleld-dependent activation en- 
ergy should be usedjASi??, Note that the obtained B{x,t) 
shown in Fig. [5] is qualitatively different from the solu- 
tion of Eq. ([2]) when the field is switched on at t = 0. 
In that case the magnetic field enters the sample in the 
form of a flux front propagating from the edges i22i22ii22 
Magnetic induction in the vicinity of the front is a linear 
function of the coordinate, B{x,t) = (An-q/^o) Vf\x — Xf\, 
with Xf(t) and Vf being the position and the velocity of 
the front. In our case, the field is switched off at i = 
and the flux escapes the sample with no front in B{x,t) 
formed even at i <C tq. Indeed, at a sufflciently large 
distance u from the edge, Eq. ^ can be linearized with 
respect to SB = Bq — B, which gives the exponential de- 
cay 6B{u,t) oc {u/2KVt)-^ ^-u^/4Kh ^ ^$oBo/47r77) 
characteristic of the linear diffusion. 

In the following we study the dynamics of the aver- 
age magnetic induction B{t) = J dxdy B{x,y,t) {A 
is the sample area) which is proportional to the mag- 
netic moment m(f) that can be measured. There are 
two regimes of the flux dynamics in the system. At 
very short times t <^ tq after switching off the fleld, 
the flux flow is localized near the edges and is unaf- 
fected by the sample size. In this case, the solution for 
a half-infinite superconductor is a good approximation, 
B{u, t) = Bof{u/ K\/t)M Here, / is a dimensionless func- 
tion to be determined from Eq. ([2]) for the half-inflnite su- 
perconductor with the boundary conditions /(O) = and 
/(oo) — 1. Using the above expression for B{u,t) and 
taking into account that it deviates signiflcantly from Bq 
in the vicinity of the edges, we find for the average induc- 
tion 1 — B[t)/ Bq cx {P/A)K\/t, where P is the perimeter 
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FIG. 4. Experimental data-'^* (o) for the magnetic moment 
m{t) in BSCO single crystal fitted to Eq. ([9]) with mo = 1.1 x 
10"^Am^ and r = 0.43 min (solid curve). This corresponds 
to »7 = 0.5Ns/m^. Inset: The inverse magnetic moment as 
a function of time. The crossover between flux flow and flux 
creep regimes is seen as a dramatic change of the slope at 
t « 7min (dashed line). The sample is a slab 0.14 x 1.37 x 
2.06 mm in size, the initial magnetic induction Bo = 35 mT, 
and T = 77K. 



of the sample. This gives the magnetic moment relax- 
ation 



m{t) = mo (^1 - y/t/f^ , 
with the time constant 

f = 9n{A/Pf7^/<PoBo 



(7) 



(8) 



where the numerical prefactor characterizes the spatial 
spread of B away from the edges. Comparison with the 
numerical solution for m{t) is shown in the inset of Fig. [3] 
for different sample geometries. We find that Eq. ([7|) is a 
good approximation of the exact m{t) in the short initial 
time interval t/f < 0.1 before the flux flow reaches the 
center of the sample. The flux flow in this time interval 
is very fast, leading to a 30% reduction of the overall 
magnetic moment. 

At times t ^ tq the flux flow extends through the 
whole sample, giving rise to the magnetization relaxation 
which depends on geometry. For the superconducting 
slab, the first-order approximation of Eqs. (HI for fc = 1 
reads B^^\t) = {%Bq/t:'^)[1 -f t/(0.75ro)]-i. Truncating 
Eqs. (Ill) at fc ^ 40, a practically exact solution is ob- 
tained. This solution can be approximated by a simple 
formula, B{t) — Bq[1 -1- t/(0.62ro)]^^, which is very close 
to the exact one for i > tq. This suggests that the exact 
solution for the magnetic moment m{t) can be accurately 
approximated by 



■m{t) = 



mo 



1 + t/i 



T = aro, 



(9) 



where a is a number which depends on geometry. Fitting 
the exact numerical solution for m{t) to Eq. ^ we find 
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FIG. 5. Experimental data^ (o) for magnetic moment re- 
laxation in NbSe2 monocrystal fitted to Eq. ((7]) (solid curve) 
with f — 1.13 X 10^ min, corresponding to = llNs/m^. The 
dashed line indicates a crossover between the flux flow regime, 
Eq. IT]), and the slow quasistatic motion before the flux creep. 
The sample has a square geometry 0.5 x 0.5 x 0.2 mm, the ini- 
tial magnetic induction Bo = 3.3 mT, and T = 4.2K. 



a = 0.620, 0.244, and 0.226 for the slab, square, and disk 
geometries, respectively (Fig. [3]). The fitting ensures the 
smallest absolute error between exact and fitted m(i) for 
1 < t/r < 3. As expected, the decay of m{t) is slower 
(that is, geometric factor a is larger) for the slab than 
for the disk, other parameters being equal. 

In what follows, we analyze available experimental 
data on m{t) and extract the characteristic time constant 
as well as the effective drag ij. Relaxation of the magnetic 
moment in BSCO single crystals is studied in Ref.[lJ. Ex- 
perimental data are shown in Fig. U] (open circles) fitted 
to Eq. dni) (sohd curve) with mg = 1.1 x 10~^Am^ and 
T = 0.43 min. The fitting is performed for the initial time 
interval before logarithmic in time, thermally activated 
fiux creep sets in. The linear time dependence of the in- 
verse magnetic moment is shown in the inset of Fig. |4l 
the crossover between fiux flow and flux creep regimes is 
seen as a dramatic change of the slope at t/r k, 16. 

Let us now extract rj. The dimensions of the sample 
used in the experiment are 0.14 x 1.37 x 2.06 mm, which 
gives Bq = A'Kmo/V = 35 mT, where V is the volume. 
Taking a = 0.620 for the slab of the width L = 0.14 mm, 
we obtain = 0.5 Ns/m^. This value for the effective 
vortex viscosity exceeds by six orders of magnitude the 
Bardeen-Stephen drag coefficient 770 ~ 10~^Ns/m^ mea- 
sured in BSCO.^- Huge enhancement of the drag indi- 
cates a strong infiuence of the pinning on the vortex dif- 
fusion. Despite the complexity of the vortex motion in 
the presence of pinning and thermal agitation, the mag- 
netization follows a simple algebraic time dependence, 
Eq. (0). 

Vortex dynamics has been studied in NbSe2 using the 
decoration technique for visualization of fiowing vortex 
lattices. Magnetization measurements have been per- 
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FIG. 6. Experimental datar^ (o) for magnetic moment relax- 
ation in YBCO polycrystal at temperatures T — 30, 41, 50, 
61, and 77 K (top to bottom), fitted to Eq. (O (solid curves) 
with f = 2.1x 10^ 5.0 X 10'', 1.4x10*, 3.0x 10^ and 7.4 x 10^ 
min, respectively. This corresponds to the drag coefficients 
ri = 284, 66, 17, 2.93, and 0.19 Ns/m^. Inset: Logarithm of 
the drag coefficient, normalized to 7;(77K) = 0.19 Ns/m^, as 
a function of the temperature. The sample has a rectangular 
geometry of 66 x 34 x 15 mm. The initial magnetic induction 
is Bo = 3.95, 3.77, 3.48, 2.80, 0.735 T, respectively. 



formed using the SQUID (superconducting quantum in- 
terference device) magnetometry. A crossover has been 
observed as a function of increasing flux density from 
a layered (smectic) flowing flux lattice in the disorder- 
dominated low-field limit to a more ordered (Bragg glass) 
lattice structure in the interaction-dominated high-field 
case. The observed time dependence of magnetization 
relaxation in the high-field limit (i?o = 3.3 mT) is shown 
in Fig. [5l The regimes indicated in Fig. [5] correspond to 
the flux flow and to the quasistatic vortex motion. The 
solid curve in Fig. [5] is the flt of m{t) to Eq. ([7]) for the 
NbSe2 sample 0.5 x 0.5 x 0.2 mm in size, which gives the 
relaxation time f = 1.13 x 10^ min and the viscous drag 
coefficient ry = 1 1 Ns/ m^ . We observe that the simple 
hydrodynamic model with an effective viscous drag force 
fits the data in the initial stages of magnetization relax- 
ation where the vortex density is large and the flux flow 
takes place. The flux flow is localized near the edges, 
as corroborated experimentally by a small reduction of 
the magnetic moment of the sample over the measure- 
ment time and, more directly, by observing the static 
vortex structure in the center of the sample.™ Large ef- 
fective rj is clearly due to hopping caused by successive 
pinning and thermally-assisted depinning of vortices, as 
evidenced by studying single- vortex dynamics in pristine 
NbSe2 monocrystals by scanning tunneling microscopyJ^ 

Magnetic moment relaxations in YBCO polycrystalii 
and monocrystaliS are shown in Figs. [6] and [T] The relax- 
ation in the polycrystalline YBCO (rectangular geome- 
try, 66 X 34 X 15mm in size) is studied at 30, 41, 50, 
61, and 77K. The initial stage of magnetization relax- 
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FIG. 7. Experimental data^* (o) for magnetic moment relax- 
ation in YBCO monocrystal at temperatures of 85 and 87 K 
(top to bottom), fitted to Eq. with f = 3.8 x 10^ min and 
f = 570 min, respectively (solid curves). The corresponding 
drag coefficients are 77 — 0.27 and 0.04 Ns/m^. The sample 
has a square geometry of 1 x 1 x 0.02 mm, and the initial 
magnetic induction is Bq =0.1 mT. 



ation can be fitted by Eq. ([7]) describing the flux flow in 
the vicinity of the edges (Fig. [SJ solid curves). This is 
in agreement with the observed small reduction of the 
overall magnetic moment during the measurement. The 
obtained effective viscosity strongly depends on tempera- 
ture, ranging between ~ 100 and 0.1 Ns/m^ as the tem- 
perature is increased from 30K to 77K, see inset of Fig. [HI 
The extracted value 7/(77K) = 0.19 Ns/m^ is consistent 
with the value 77 — 0.12 Ns/m^ measured independently 
at the same temperature by studying the spatiotempo- 
ral change of the magnetization profile in a bulk YBCO 
sample in the flux-flow regime. Taking 77 oc ^ 1^"^ and 
neglecting the temperature dependence of the effective 
activation energy V as well as of the prefactor, we find 
U « 360 K in accordance with the previous results<^ 

Magnetic moment relaxation in small YBCO 
monocrystal (1 x 1 x 0.02 mm) is shown in Fig. [71— 
The data can be fitted with the effective viscous drag 
coefficients 77 = 0.27 and 0.04 Ns/m^ at temperatures 
of 85K and 87K, respectively. The decrease of 77 in 
such a narrow temperature range may be due to the 
proximity of the critical temperature (Tc ~ 88K) where 
fluctuations are pronounced. In addition, the sharp 
change in the relaxation rate observed at 87K and 
t w 30 min suggests that the flux flow is inhomogeneous 
and made of large domains which, upon depinning, 
abruptly increase the magnetic moment relaxation rate. 

In conclusion, we have studied vortex dynamics in 
type-II superconductors in the initial time interval be- 
fore the flux creep sets in. We have used a simple phe- 
nomenological (hydrodynamic) model of nonlinear diffu- 
sion of massless vortices where pinning of the flux lines by 
material inhomogeneities, interaction with other vortices 
and the surface, and the Bardeen-Stephen dissipation in 
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the vortex core are described by an effective viscous drag 
coefficient, rj. After switching off the external magnetic 
field, the vortex dynamics exhibits two distinct regimes 
before the creep sets in with logarithmic in time decay 
of remanent magnetization. In the beginning, the flux 
flow is localized near the edges and is independent of the; 
sample size. At later times, this regime is followed by 
a slower flux flow involving the bulk of the sample. We 
find that magnetic moment relaxation in these regimes 
can be accurately approximated by m{t) = mo(l— -^/t/f) 
for f ^ f and m{t) = mo(l + t/T)~^ for t > t, where 
geometry-dependent f and t arc proportional to 77. 

We have analyzed available experimental data on 
early stages of magnetization relaxation after the mag- 
netic field is instantaneously removed. We obtained 
quantitative agreement for both conventional and high- 
temperature superconductors, albeit with rj exceeding 
the Bardeen-Stephen value % by many orders of magni- 



tude. Huge enhancement of rj with respect to 770, as well 
as its exponential temperature dependence, indicates a 
strong influence of pinning and thermally assisted depin- 
ning of vortices on flux diffusion. We argue that early 
stages of magnetization relaxation can be modeled as a 
flux flow with an effective drag coefficient. This allows 
for a simple experimental determination of the bulk vor- 
tex viscosity, which cannot be accessed by the surface 
impedance measurements. 
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